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Abstract 

We discuss the possibility that the PeV neutrinos recently observed by 
IceCube are produced by the interactions of extragalactic cosmic rays during 
their propagation through the radiation backgrounds. We show that the fluxes 
resulting from the decays of neutrons produced in the interactions of cosmic 
ray protons with the CMB background are suppressed (E^d$> v /dE < 10 -10 
GeV/cm 2 s sr), with those resulting from the decays of pions produced in the 
interactions with the UV/optical/IR backgrounds being the dominant ones 
at PeV energies. The anti- neutrino fluxes produced by the decay of neu- 
trons resulting from the photodisintegration of heavy nuclei with CMB pho- 
tons are also shown to be quite suppressed (E%d$> u /dE < 1CP 11 GeV/cm 2 s 
sr), while those produced by photo-pion processes with UV/optical/IR back- 
grounds may be larger, although they are not expected to be above those 
achievable in the pure proton case. Scenarios with mixed composition and 
low cutoff rigidities can lead to PeV neutrino fluxes enhanced with respect 
to those in the pure Fe scenarios. We also discuss the possible impact of the 
Glashow resonance for the detection of these scenarios, showing that it plays 
a moderate role. 



1 Introduction 

The IceCube Collaboration recently reported [1] the observation of two neutrinos 
with energies in the range 1-10 PeV, which are the highest energy neutrino events 
observed up to now. Possible sources for their origin include atmospheric neutrinos 
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(in particular from the prompt charm production), astrophysical neutrinos produced 
by photo-pion or pp interactions of accelerated protons with ambient radiation or 
matter in sources such as gamma ray bursts, active galactic nuclei or starburst 
galaxies, or cosmogenic neutrinos produced during the propagation of ultra-high 
energy cosmic rays (UHECR) through the extragalactic radiation backgrounds. 

Although two events are certainly not enough to determine the neutrino flux, 
these observations suggest a flux level of E%d$ v /dE ~ 1CT 8 GeV/cm 2 s sr in the 
1-10 PeV range [1]. Note that the canonical expectations for the rapidly falling 
atmospheric background neutrino flux are about one order of magnitude below this 
level, while astrophysical source scenarios may produce the required flux. It is our 
purpose here to discuss some generic upper bounds on the cosmogenic neutrino 
fluxes of various origins to identify the potentially dominant contributions at PeV 
energies. 

Another relevant feature of the observations is that the two neutrinos are in the 
cascade mode (i.e. not involving muon or tau tracks), so that they are due to either 
a NC interaction of any of the three neutrino (or antineutrino) flavors or to charged 
current interactions of electron neutrinos (or antineutrinos). One interesting aspect 
of this is that the P e channel presents, besides the interaction with nucleons, a 
resonant interaction with the electrons at an energy of Eqr = My V /2m e = 6.3 PeV 
(the Glashow resonance [2]), and we will also discuss the possible impact of this 
channel for the IceCube detection. 

2 Neutrino fluxes from UHECR protons 

Extragalactic UHECR protons with energies above 6 x 10 19 eV get attenuated 
when propagating through the cosmic microwave background (CMB) mostly by 
pion production processes, which lead to the well known expected GZK suppres- 
sion [3]. These losses reach a maximum strength at the A(1232) resonance, i.e. 
when m\ ~ m 2 + 2E p E 1 {l — cos 9), with 9 the angle between the p and 7 mo- 
menta in the lab frame. This corresponds to a proton energy (adopting 9 = ir) of 
E P (A) ~ 1.6 x 10 20 eV/(£ 7 /l(T 3 eV), with the average CMB photon energy being 
~ 0.7 x 10 -3 eV at present (i.e. for redshift z = 0). The process pq — > n + n leads 
to the production of cosmogenic neutrinos [4] both through the pion decay chain 
7r + — > —> e + u e u fl u^ l and through the neutron decay n — > peu e [5]. Since the 

pions typically carry about 1/5 of the proton energy, each neutrino from the pion 
decay has on average an energy of about E v ~ E p /20. On the other hand, the v e 
from the neutron decay has a typical energy of E v ~ 4 x 10~ 4 E n ~ 3 x lO^ 4 ^. 

When these interactions happen at high redshift, since the CMB temperature 
scales as (1 + z) the proton energies for which the photopion production start to be 
efficient are E p ~ E GZK /(1 + z), where for definiteness we adopt E GZK = 10 20 eV 
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as the typical proton energy for pion production at z — 0. Note also that the 
CMB photon density increases as (1 + z) 3 , making the opacity of the universe to 
UHE protons correspondingly higher. The energies of the neutrinos produced at 
high redshifts get further reduced by the adiabatic losses as they propagate to us, 
leading to E v ~ Egzk/{2®0- + z ) 2 ) f° r the neutrinos originating from pion decays, 



and E v 



3 x 10" 



Egzk/ (1 



for those from neutron decays. The neutrino 



production turns out to be sizeable up to redshifts of 3-5 [5] , depending on the actual 
source redshift evolution, and for instance considering a typical neutrino production 
redshift z ~ 1.2 one gets peaks on the neutrino spectrum resulting from interactions 
with the CMB at energies ~ 10 18 eV (from 7r decays) and ~ 6 x 10 15 eV (from 
n decays). These peaks are anyway quite wide, because the A resonance is wide 
(and other pion production channels contribute as well), because the CMB photons 
have a wide thermal spectrum and also because different redshifts contribute to the 
neutrino production. 

A useful relation can be obtained between the two neutrino fluxes just considered, 
since in the charged pion producing interactions the same number of low energy v e 
and higher energy u e , or are produced. Hence, denoting by § u the resulting 
neutrino diffuse fluxes, one has, for the fluxes produced in interactions with the 
CMB alone, that 
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Using that the EeV neutrinos are actually dominated by those produced in in- 
teractions with the CMB (see below), and ignoring for the time being the effects of 
the neutrino oscillations on the P e , we then get 
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The all flavor diffuse neutrino flux has been constrained at EeV energies by 
the unsuccessful searches by IceCube [6] and Auger [7], which imply the approxi- 
mate bound E^d&v/dE < 3 x 10~ 7 GeV/cm 2 s sr at EeV energies. Moreover, in 
proton scenarios a stronger bound has been obtained indirectly from the so-called 
cascade decays [15]. This bound is derived from the requirement that the 7r° de- 
cay gammas (also produced in the photo-pion processes) and the e + e~ pairs do 
not produce too large amounts of GeV-TeV photons when they cascade down to 
low energies as they interact with the CMB and IR radiation backgrounds. The 
allowed amount of low energy photons is bounded by the diffuse photon back- 
ground measured by the Fermi LAT experiment [16]. The cascade limit which 
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results for the all flavor cosmogenic neutrino flux at EeV energies is at the level of 
E^dfyv/dE < 5 x 10~ 8 GeV/cm 2 s sr [17, 18]. This, combined with eq. (2), implies 
that the v e flux at PeV energies produced from interactions with CMB photons 
should satisfy E^dQpJdE < 10~ 10 GeV/cm 2 s sr. This upper bound is about two 
orders of magnitude below the flux level suggested by the two PeV neutrinos ob- 
served by IceCube, and hence can hardly be responsible for those events (this is at 
variance with the interpretation suggested in [19]). 

Cosmic ray protons with energies below 6 x 10 19 eV (at z — 0) loose energy mainly 
by e + e~ pair creation interactions with CMB photons, without producing neutrinos. 
However, below ~ 10 18 eV the dominant attenuation process becomes the photopion 
production in interactions with UV, optical and IR radiation backgrounds, which 
being more energetic than CMB photons lead to a reduced proton energy to produce 
the A resonance (E p (A) ~ 1.6 x 10 17 eV/(£' 7 /eV)). In this case, the neutrinos 
produced by the pion decays have energies ~ 8 x 10 15 eV/[(l + z)(E 7 /eV)], and 
hence will contribute in the PeV range, while those from the corresponding neutron 
decays will be at much lower energies (< 10 14 eV), and hence buried below the 
atmospheric neutrino background. The PeV neutrino fluxes from UHECR proton 
scenarios are indeed dominated by those from pion decays produced in interactions 
with UV/optical/IR photons. Numerical simulations show that all flavor neutrino 
flux levels of E^d^ u /dE ~ few x 10 _9 GeV/cm 2 s sr can be achieved for E v ~ PeV 
(see e.g. [9, 10]), with the precise values depending on the assumed source evolution, 
on the radiation background adopted and on the shape of the extragalactic proton 
spectrum at the source. 

To illustrate this, we show in figure 1 the neutrino fluxes resulting in pure proton 
scenarios, adopting a source power spectrum a = 2.4, maximum energy of 200 EeV 
and minimum energy of 2 x 10 16 eV. We also adopt a source redshift evolution (for 
the density times CR emissivity) following the gamma ray bursts one (correspond- 
ing to the SFR6 model derived in ref. [20], here referred to as GRB2). For the 
UV/optical/IR radiation background we consider (in all the figures) the one follow- 
ing the 'best fit model' in ref. [21] (including also its redshift evolution). We obtain 
the CR spectra as well as the diffuse neutrino fluxes using the simulation code CR- 
Propa [23]. Besides the total fluxes, we show separately the contributions resulting 
from the interactions with CMB photons and that of antineutrinos from n-decays 
alone. The relation obtained in eq. (2) can be seen to hold by comparing the heights 
of the two peaks (the EeV all flavor one and the PeV one from n-decays). 

The resulting CR spectrum is normalized at 10 EeV to that measured by the 
Auger Collaboration [24] (which has a 22% systematic uncertainty on the energy 
scale, not shown in the plot). The overall shape of the spectrum is found to be 
in reasonable agreement with the measured one above the ankle region. We also 
display the spectrum measured by the Hires Collaboration [25], and the bounds on 
the all flavor neutrino fluxes obtained by IceCube [6], Auger [7] and Anita [26]. 
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proton sources, oc=2.4, E max =200 EeV, GRB2 




Figure 1: Proton 'dip' scenario with source spectral index a = 2.4 and E max = 
200 EeV. Indicated are the propagated proton spectrum and the resulting all flavor 
neutrino fluxes (obtained with CRPropa). We also show separately the neutrino 
backgrounds due to interactions with CMB alone as well as those resulting from n 
decays. The CR flux measured by Auger [24] and Hires [25] and the neutrino limits 
from IceCube [6], Auger [7] and Anita [26] are displayed. We also indicate the energy 
range and approximate flux level suggested by the two observed IceCube events. 

The PeV neutrino flavors produced in these scenarios, arising mostly from pion 
decays, are in the ratio {y e : : v T ) = (1 : 1 : 0) and {p e : : v T ) = (0 : 1 : 0). 
These ratios will then be affected by the incoherent flavor oscillations that take place 
from their production until their detection on Earth. For instance, if one adopts 
the tri-bi maximal (TBM) neutrino mixing pattern, these fluxes get transformed 
into (0.78:0.61:0.61) and (0.22:0.39:0.39) upon arrival on Earth [13]. Departures 
from the TBM mixing predictions, as required by the recent measurement of a non- 
vanishing # 13 [14], will induce small shifts on the above mentioned flavor ratios. Note 
that in the case of v e production in neutron decays, the initial flavor ratios of the 
antineutrinos (1:0:0) would change into (0.56:0.22:0.22) by oscillations in the TBM 
scheme. 

In figure 2 we show the predicted CR spectra, the all flavor neutrino fluxes 
as well as the cascade photons for three different source evolution models. These 
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Figure 2: Proton scenario with E max = 200 EeV for different source evolution models 
(SFR1, GRB2 and FRII). The source spectral index is a = 2.4 for the SFR1 and 
GRB2 models, while a = 2.2 for the FRII model. Indicated are the propagated 
proton spectrum, the resulting (all flavor) neutrino and the photon fluxes. The 
photon background measured by Fermi-LAT [16] is indicated, besides the CR spectra 
and v bounds included in fig. 1. 

models follow the star formation rate (SFR1 in [20]), the GRB2 model considered 
in fig. 1 and the AGN model evolution from ref. [22] (FRII). In order to obtain a 
reasonable fit to the high energy cosmic ray spectrum, we adopt a = 2.4 for the 
SFR1 and GRB2 models, while a = 2.2 for the FRII one. A significant dependence 
of the predictions on these parameters is observed. We see that the gamma ray 
fluxes obtained in the scenario with the strongest neutrino fluxes are on the verge 
of conflicting with the diffuse gamma ray flux observed by Fermi [16]. 

One further issue is that the height of the PeV neutrino peak (resulting mostly 
from interactions of ~ 10 17 eV protons with UV/optical/IR photons) and that of 
the EeV peak (resulting from interactions of ~ 10 20 eV protons with CMB photons) 
are in principle related, depending on the underlying proton spectral shape. Hence, 
the bounds on the EeV diffuse neutrino fluxes (direct or from cascade decays) may 
also constrain the allowed maximum height of the PeV neutrino peak. We note 
however that if the proton spectrum above a few EeV (or even a few tens of EeV) 
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were suppressed, for instance due to the existence of a low maximum rigidity in 
the acceleration process at the sources, the EeV neutrino peak would be suppressed 
but not the PeV peak, and hence in these cases no constraint from the associated 
cascade decays would result. 

Of interest in this respect is the so-called disappointing model [27], which actually 
has an enhanced proton flux at EeV energies, with a cutoff at few EeV, and in this 
scenario the cosmic rays above the ankle at ~ 4 EeV are nuclei with increasingly 
heavier masses. This model would predict an enhanced neutrino flux at PeV energies 
with no sizeable neutrino peak at EeV energies (see next Section). 

We also note that the PeV neutrino flux is sensitive to the amount of UV photons, 
and hence scenarios with an enhanced UV background, as considered in [9, 10], can 
lead to some enhancement in PeV neutrino fluxes. 



3 Neutrino fluxes from UHECR nuclei 

Scenarios in which heavier nuclei make a significant contribution to the UHECR flux 
are qualitatively different. Here the photopion production off the CMB photons only 
occurs for energies above ~ AEqzk, where A is the mass number of the nucleus. 
On the other hand, photodisintegration processes play an important role at lower 
energies. Photodisintegrations are dominated by the giant dipole resonance (GDR) 
which, in the nucleus rest frame, has a threshold for photon energies between a few 
MeV and 10 MeV (depending on the nucleus), and peaking at about 20 MeV. The 
photon energy in the CR rest frame can be expressed as 

*• - * - - v (!) («) (^H ife < 3 > 

where (3 = \Jl — (m p A/E) 2 is the boost factor and 6 is the angle between the CR 
and photon momenta in the lab frame. While they photodisintegrate, nuclei emit 
nucleons without changing their Lorentz factor, hence just loosing energy because 
the mass of the leading fragment gets reduced. In this process secondary nucleons 
are emitted with energies E/A for energies E > Eqdr — 10 20 eVA/56. If the photo- 
disintegrations take place instead at non-zero redshift, the secondary nucleons final 
energy will be E/A{1 + z), where the contributions arise from E > E GDR /{1 + z). 
About half of these nucleons will be emitted as neutrons, which then decay to pro- 
duce V e with energies E v > (4 x 10~ 4 )(10 20 eV/56)/(l + z) 2 ~ few 10 14 eV (assuming 
interactions with the CMB background alone). Note that for interactions with the 
CMB, only those neutrinos produced by energetic nuclei at or above the GDR will 
end up with energies above one PeV. The nuclei of smaller energies may also photo- 
disintegrate by interacting with the UV/optical/IR radiation backgrounds, leading 
to a flux of v e of low energy (< 10 14 eV), and hence are not relevant here. 
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One may obtain a conservative upper bound on the flux of the v e produced at 
PeV energies by noting that the secondary nucleons produced by photodisintegra- 
tion with CMB photons end up, regardless of the primary nucleus mass, piled up 
around a few EeV, and only about half of these were produced as neutrons. Hence, 
by requiring that the secondary nucleons flux at the relevant energy be below the 
actually measured CR flux one ends up with 
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where we used that the measured CR flux at 2.5 EeV is about d^cn/dE ~ 7 x 
10~ 18 /(EeV cm 2 s sr) [24]. These bounds may be slightly modified if one considers 
that the secondary nucleon energies could be degraded by pair production processes, 
but at few EeV energies this process does not have a large impact. Hence, the v e flux 
produced by decays of neutrons resulting from the photodisintegration of UHECR 
nuclei are quite suppressed at PeV energies. 

On the other hand, neutrinos with few PeV energies may be produced by the 
interactions of extragalactic nuclei with UV/optical/IR photons by photopion pro- 
duction. This would require that the energy per nucleon be about ~ 20E U , and 
hence the nuclei should have energies of about 20 AE V . However, since the UHECR 
source spectrum is expected to be steeper than E~ 2 , the photopion contribution 
from scenarios dominated by nuclei will not be larger than that arising in the pro- 
ton dominated scenarios discussed in the previous section (for instance, if all CR 
had a given mass number A, the number of nucleons per logarithmic energy interval 
at energy E is related to the corresponding number in a pure proton scenario by the 
factor A 2 ~ a , where a is the assumed power spectrum slope). Hence, nuclei scenarios 
give rise to PeV neutrinos by photopion production of UV/optical/IR photons but 
at a level which is not expected to be larger than that achievable in proton scenarios. 

To illustrate the predictions from a heavy composition scenario we show in fig. 3 
a case corresponding to Fe only sources with a source spectrum a = 2.0 and a 
maximum Fe energy of 5200 EeV (corresponding to rigidities R = E/Z< 200 EV, 
with Z the charge of the nucleus), following the GRB2 source redshift evolution. We 
also show separately the neutrino fluxes arising from the neutron decays, where the 
bound from eq. (4) can be seen to hold. The main contribution to the PeV neutrino 
fluxes arises from the interactions with UV/optical/IR radiation backgrounds. 

Figure 4 shows instead the results obtained in a scenario having a mixture of p 
and Fe and a low energy cutoff (E max = 5Z EeV), inspired in the 'disappointing' 
model [27]. The relative source abundances considered are n v jnp e = 10 at a given 
energy (below the proton cutoff), with power spectrum a = 2.0 and for the GRB2 
source evolution model. We see that the enhanced proton contribution below the 
ankle helps to reach a larger flux of PeV neutrinos than in the pure Fe case. In 
this mixed composition scenario the CR spectrum in the ankle region is similar to 
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Fe sources, oc=2.0, E =5200 EeV, GRB2 
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Figure 3: Extragalactic Fe scenario with source spectral index a = 2.0 and E max = 
5200 EeV. Indicated are the propagated CR spectrum and the resulting (all flavor) 
neutrino fluxes, as well as the neutrino background due to n-decays alone. 

the measured one, but it does not fit well the highest energies. This may however 
depend on the precise distribution of nearby sources and on the shape of the source 
cutoff adopted. 

4 The impact of the Glashow resonance 

An interesting aspect of PeV neutrinos is that the v e may be detected through the 
interactions with electrons [2, 28, 29], which gets resonantly enhanced at the W pole, 
corresponding to an energy Eqr = My V /2m e = 6.3 PeV. The W decays with ~ 10% 
branching into each lepton flavor and the remaining 70% into hadronic channels. 
The electron and hadronic channels (i.e. about 80% of the decays) would give rise 
to a cascade signature, while the muon and tau channels give rise to tracks (or 
lollipops, since the tau has a track length of about 300 m at 6 PeV, decaying then 
to produce a cascade). The resonant cross section is (considering e.g. the muonic 
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Figure 4: Mixed composition (p-Fe) scenario with source spectral index a = 2.0 
and E max = 5 Z EeV. Indicated are the propagated CR spectrum, the resulting (all 
flavor) neutrino fluxes and the separate contributions from p and Fe primaries. 



channel) 



GR(- . - \ G 2 F Mw E V /E GR 
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The total, i.e. into all fermion final states, cross section reaches a peak value of 
5 x 10~ 31 cm 2 , having an associated width A = (Tw/Mw)Egr — 0.17 PeV, which 
is quite narrow. 

On the other hand, PeV neutrinos of all flavors interact with nucleons through 
charged current interactions with a cross section [30] 

<r cc (uiN -)• i t X) ~ l.Q(E/E GR )°- m3 10~ 33 cm 2 , (6) 

and through neutral currents with 

a N C{y . N ^ V . X) „ 6 . 8(jE ; /jE ; Gi?) 0.363 10 -34 (?) 

with the cross sections for the antineutrinos having similar values. 

To compare the relative importance of the resonant cross section on the rates 
one has to take into account the relative abundances of the different neutrino and 
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antineutrino flavors, the fact that the ratio of the electron and nucleon densities 
in water is n e /njv = 5/9 and that the expected fluxes span a wide energy range. 
Since, as shown in the previous sections, the dominant fluxes of PeV neutrinos 
are those resulting from pion decays (produced in interactions with UV/optical/IR 
photons), the flavor ratios on Earth are £ ~ (0.78 : 0.61 : 0.61) for neutrinos and 
£ ~ (0.22 : 0.39 : 0.39) for antineutrinos (adopting TBM mixing for simplicity). 
Focusing on the cascade producing channels only, the relative contribution of the 
GR channel is, assuming for definiteness that in the energy interval [E min , E max ] 
around E GR the flux spectral shape can be approximated as oc E~ 2 , 



n e Jp max dE E-%8a GR (u e e ->■ z/ M u) 

n N j§ ::: dE e-* (Ei& + iM NC + (Ce + q° cc 



From this expression one gets for instance that for a very narrow interval around 
E GR with width A = 0.17 PeV the ratio is R(E GR - A, E GR + A) = 10.7, while for 
an interval with 2.5 PeV width one gets #(3.8 PeV,8.8 PeV) = 1.23. Note that in 
the v e CC channel all the neutrino energy is transferred to the cascade, while for 
the NC ones the average inelasticity of the interaction is (y) ~ 0.26, and hence the 
cascades produced will have a lower energy, but for the range of neutrino energies 
considered they will still be above PeV energies. Note that the ratio of resonant 
to non-resonant rates further decreases for wider energy intervals or if one were to 
include the muon and tau channels. 

We see that although the total peak cross section at E GR is ~ 350 times the value 
of a cc , taking all the relevant factors into account the final enhancement which can 
be obtained with this process for a wide neutrino spectrum around E GR is modest 
(of order unity) 1 . 



5 Conclusions 

We have discussed the impact of the different production mechanisms of cosmogenic 
neutrinos on the fluxes at PeV energies. The main results are that: 

- The neutrinos resulting from the decays of neutrons produced in photo-pion 
interactions off the CMB are at a level Eld<$> u /dE < 10" 10 GeV/cm 2 s sr at PeV 
energies, and hence make a negligible contribution. 

- The neutrinos resulting from the decays of neutrons produced by photo-disintegration 
of heavy nuclei are at a level Eld$ v /dE < 10~ n GeV/cm 2 s sr, and are hence also 
negligible. 

- The neutrinos resulting from pion decays produced in photo-pion interactions 
with UV/optical/IR photons are the dominant ones at PeV energies, and can reach 

1 While this work was being finished, a preprint appeared discussing possible signatures of the 
Glashow resonance and non-standard scenarios to account for the two IceCube events [31]. 
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the level of E^&<& V /&E ~ 10~ 9 GeV/cm 2 s sr in proton scenarios, and somewhat 
lower in heavy nuclei scenarios. 

- The PeV neutrino fluxes may be enhanced in scenarios with strong source 
redshift evolution or steeper source spectra, and also if the UV background were 
strongly enhanced at high redshift. 

- The overall contribution to the observable rates in the 1-10 PeV range due 
to the Glashow resonance is not large for the cosmogenic neutrinos. For instance, 
the v e e resonant interaction leads to a contribution to the rates of cascade events 
comparable to the non-resonant one resulting from the neutrinos within 2.5 PeV of 
Egr- 

When taking into account the CR primary spectrum as well as the secondary 
gamma ray cascade flux it turns out to be difficult to interpret the two IceCube 
events in terms of a cosmogenic neutrino flux, unless the two events were an upward 
fluctuation of the expected neutrino rates. Hence, the interpretation in terms of 
neutrino production at the sources seems likely, and these need not be necessarily 
at high redshifts. 
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